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ABSTRACT 

We calculate detailed NLTE synthetic spectra of a Pulsating Reverse Detona- 
tion (PRD) model, a novel explosion mechanism for Type la supernovae. While 
the hydro models are calculated in 3-D, the spectra use an angle averaged hydro 
model and thus some of the 3-D details are lost, but the overall average should 
be a good representation of the average observed spectra. We study the model at 
3 epochs: maximum light, seven days prior to maximum light, and 5 days after 
maximum light. At maximum the defining Si II feature is prominent, but there 
is also a prominent C II feature, not usually observed in normal SNe la near 
maximum. We compare to the early spectrum of SN 2006D which did show a 
prominent C II feature, but the fit to the observations is not compelling. Finally 
we compare to the post-maximum UV+optical spectrum of SN 1992A. With the 
broad spectral coverage it is clear that the iron-peak elements on the outside of 
the model push too much flux to the red and thus the particular PRD realiza- 
tions studied would be intrinsically far redder than observed SNe la. We briefly 
discuss variations that could improve future PRD models. 
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Introduction 



Understanding the nature of the explosion mechanism of Type la supernovae (SNe la) 
is important for the use of SNe la as precision cosmological probes as well as for stellar 
evolution and nucleosynthesis. Determining something about the nature of the explosion 
will almost certainly shed ligh t on the phys i cal explanation o f the e mpirical light curve 



shape-luminosity rela t ionship Jphillipsl 



Goldhaber et al.1 12OOII : iKasen fc Wooslev 



993 



Riess et all Il995l . Il996l : IPhilhps et all Il999 



20071 ) as well as further physical parameters as- 



sociated with the observed diversity in the SN la population. Understanding the physical 
nature of the diversity is crucial in being able to estimate systematic errors that may arise 
in comparing the characteristics of the nearby sample with the cosmological one. 

While there is widespread agreement that the SN la progenitor consists of a C+0 white 
dwarf that accretes material from a companion and ignites the C+0 fuel when the Chan- 
drasekhar mass is reached, the exact mechanism of explosion is not known. The explosion 
could take the form of a super-sonic detonation (shock- wave), or a sub-sonic deflagration. If a 
unexpanded Chandrasekha r mass WD detonat es, the resulting products would be almost en- 
tirely ^^Ni and a-particles ( lArnett et al.lll97ll ). contrary to the signature intermediate mass 
elements of silicon, sulfur, and calcium that define the SNe la class. For many years the stan- 



dard SN la model has been the 1-D, parameterized, deflagration model W7 (INomoto et al. 



19841 ). In that model the speed of the flame was adjusted in order to reproduce the approxi- 
mately correct amount of ^^Ni, as well as intermediate mass elements. Nevertheless, the cen- 
tral regions of the model undergo significant electron ca pture and far too much neutron-rich 
material is produced, polluting the interstellar medium. iBrachwitz et al.l (2000|) showed that 
with pf-shell mo del calculations the electron capture rat es are reduced over those of standard 
sd-shell rates of Fuller et al.l (jl980l : Il982bl : Il982al : Il985l ). we are unaware of hydrodynamical 
calculations that show that the problem is completely solved. In order to avo id some of the 



problems of W7, the delayed-detq nation model has gained significant favor (IGamezo et al. 



2OO5I . l2004j : iKhokhlovlEJ . Il991al lb[). In this scenario, the ignition begins as a deflagration, 
allowing the star to pre-expand, before the deflagration transforms into a detonation at lower 
densities, which produces both ^^Ni and intermediate mass elements. While this process is 
known to occur in confined terrestrial situations, it is not known whether it can take place 
in the unconfined white dwarf. 



The last decade has seen the beginning of full 3-D hydro calculations 
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Ropke et al.ii2006i: iRopke et al.ll2006: Ropke & Hillebrandt 2005a b: I 


■lopke 


2005). Much 



of the work has focused on pure deflagration models as well as developing the necessary 
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sub-grid techniques to handle the essentially infinite djTiamic range (due to turbulence) or 
studying the effects of differing numbers of ignition points. The initial ignition conditions 
may well be unki iowable since abou t one hundred years before ignition the white dwarf be- 
comes turbulent (iRopke et al.l 120071 ). Recently, there has been general agreement that pure 
deflagration models leave behind significant amounts of unburnt carbon and oxygen and that 



the kine tic energy of the exp 
SNe la JCamezo et all booi 



(iKhokhlov 



1993 



1991a 



b 



osion is too low to reproduce the observed spectra of normal 



Mazzali et a. 



Gamezo et al. 



20071). The alternative is a delayed detonation 



20051 : iMazzali et al.l 120071 ) where the defiagration 



accelerates into a detonation or some other process that turns the defiagration into a det- 



onation. Two alternatives that act in a similar manner to t 



are the "Gravitationally Confined D etonation" (jPlewa et a.1 



le delayed-detonation sce nario 



the "Pulsating Reverse Detonation" ( Bravo fc Garcia-Senzl 2006 ) 



20M; [Jordan et all 120071 ) and 



We present detailed NLTE synthetic spectral calculations of an angle averaged ver- 
sion of PRD5.5 a pulsating-reverse-detonation calculation (E. Bravo & D. Garcia-Senz, in 
preparation). The detailed chemical composition was determined from the hydro model 
with a post-processing code that included a maximum of 725 isotopes. The evolution of the 
electron fraction was followed along with the hydrodynamic calculation with the aid of 
tables of electron captur e rates in NSE matter. The procedure was identical to that used in 
Garcia-Senz et al.l (119991 ). 



2. Calculations 



The calculation s were performed using the multi-purpose stellar atmospheres program 



PHOENIX version 15 (IHauschildt fc Baronlll999l : iBaron fc Hauschildtlll998l : iHauschildt et al. 



1997al lbl. E996l ). Version 15 incorporates many changes over previous versions used for super- 
nova modeling (IBaron et al.ll2007l . l2006l ) including many more species in the equation of state 
(83 versus 40), twice as many atomic lines, and many more species treated in full NLTE. 
PHOENIX solves the radiative transfer equation along characteristic rays in spherical symme- 
try including all special relativistic effects. The non-LTE (NLTE) rate equations for many 
ionization states are solved including the effects of ionization due to non-thermal electrons 
from the 7-rays produced by the radiative decay of ^^Ni, which is produced in the supernova 
explosion. The atoms and ions calculated in NLTE are: H I, He I-II, C I-III, O I-III, Ne I, 
Na I-II, Mg I-III, Si I-III, S I-III, Ca II, Ti II, Fe I-III, Ni I-III, and Co I-III. These are all 
the elements whose features make important contributions to the observed spectral features 
in SNe la. 



Each model atom includes primary NLTE transitions, which are used to calculate the 
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level populations and opacity, and weaker secondary LTE transitions which are included in 
the opacity and im plicitly affect the rate equa tions via their effect on the solution to the 
transport equation (IHauschildt fc BaronI Il999l ). In addition to the NLTE transitions, all 
other LTE line opacities for atomic species not treated in NLTE are treated wit h the equiva- 
lent t wo-level atom source function, using a thermalization parameter, a = 0.10 (IBaron et al. 
19961 ). The atmospheres are iterated to energy balance in the co-moving frame; while we 
neglect the exphcit effects of time dependence in the radiation transport equation, we do 
implicitly include these effects, via explicitly including the rate of gamma-ray deposition in 
the generalized equation of radiative equilibrium and in the rate equations for the NLTE 
populations. 

The outer boundary condition is the total bolometric luminosity in the observer's frame. 
The inner boundary condition is that the flux at the innermost zone [v = 1500 km s~^) is 
given by the diffusion equation. We treat the 7-ray deposition by solving the full radia- 
tive transfer for a grey opacity. Detailed fitting of the observed spectra determines all the 
parameters. 

19933. 



PHOENIX has b een well tested on SNe la (iNugent et al 



Baron et al.l |2006| ) and particularly on SN 1994D (ILentz et al 



199 



I2OOI 



Lentz et al 



Baron et al 



2001 



20061). 



PHOENIX is subjected to rigorous regression testing with each new version. 



3. The Model 



In the Pulsating Reverse Detonation model (IBravo fc Garcia-Senzl l2006l ) the deflagra- 
tion starts in a number of rising plumes that are already somewhat close to the surface. 
These plumes burn without much expansion of the white dwarf and hence burn to iron peak 
elements. The energy released by this initial deflagration is not enough to unbind the star 
and the outer parts fall back producing an accretion shock. The unburnt core settles back 
into equilibrium and the accretion shock detonates the outer parts which are at low enough 
density that they burn to intermediate mass elements. Thus the composition structure in 
this model is not completely stratified and the outer layers are polluted by iron and radioac- 
tive nickel. The precise amount of polluting Fe-group elements is sensitive to the outcome of 
the defiagrative phase of the explosion, more efficient subsonic burning reduces the amount 
of stable Fe-group pollution. The mass fraction of radioactive material close to the surface 
of the ejecta depends as well on the amount of electron captures that, in turn, are a sensitive 
function of the central density at ignition. 



We present results from model PRD5.5 (E. Bravo & D. Garcia-Senz, in preparation). 
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The model started from five sparks, incinerating 0.14 Mq during the initial subsonic combus- 
tion phase, and later on a detonation processed most of the remaining carbon and oxygen. 
Finally, the explosion produced 1.3 x 10^^ ergs of kinetic energy, and ejected 0.86 Mq of 
^^Ni, 0.15 Mq of intermediate-mass elements, and 0.18 of carbon and oxygen. In model 
PRD5.5 the ejected mass of stable iron and nickel was O.IM©, most of which was moving at 
velocities in excess of 10,000 km s~^. The final averaged di stribution of elements in vel ocity 
space is quite similar to that of model PRD6 (see Fig. 3 in lBravo fc Garcfa-Seng|2006l ). 



The 3-D composition of model PRD5.5 is quite clumpy. The angle averaged version 
is rather uniform in the outer layers above ~ 8000 km s~^, where the supernova spectrum 
mostly forms before a few days past maximum light (Fig. [1]). While the outer layers are 
dominated by ^^Ni, there is significant silicon, sulfur, and oxygen as required for SN la 
spectra. It is a quantitative question that only a code such as PHOENIX can determine 
whether or not the outer ^^Ni conflicts with observations. There is also significant carbon 
in the outer laye r s and C+0 dominates the innermost region. Syn thetic spectru m modeling 
by iKozma et al.l (120051 ) of the deflagration model c3_3d_256_10s (iRopkd l2005l ) found that 
abundant C+0 matter in the interior was inconsistent with late-time nebular spectra. Model 
c3_3d_256_10s had more interior C+0 than PRD5.5, therefore future nebular modeling of 
PRD5.5 is also important. 



Results 



We first examine the maximum light spectrum of SN 1994D, which we have used in 
previous pape rs to study the detailed spectral beh avior since it is nearby and extremely 
well-observed ( iBranch et al.ll2005l : iLentz et al.ll200ll ). Figure [2] shows the results of our best- 
fit model (we use the outer boundary condition of the total bolometric luminosity in the 
observer's frame to find the best fit). The overall fit is not all that bad, the feature at 
6150 A attributed to Si II A6355 is reasonably well fit, however the 'W" feature attributed 
to S II is totally washed out. Additionally just to the red of the 6150 feature is another 
feature almost c ertainly due to C II A6580. This feature has in fact been observed in early 
SNe la spectra (IBranch et al.ll2003l : iThomas et al.l 120071 ). While the 6150 A feature is the 
usual defining feature for SN la, the S II "W" is even a more reliable diagnostic. SNe Ib/c 
spectra do show features that have been identified as Si II, but only SNe la show the S II 
feature. Since silicon and sulfur are made nearly simultaneously one might expect the lines 
to appear together, however the Si II line is much stronger than the S II lines, thus it requires 
more sulfur to obtain the S II 'W" feature. Empirically S II seems to be the best classifier 
of SNe la. Even the extremely odd SN la 2002ic which shows strong narrow Balmer lines 
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Fig. 1. — The angle averaged composition of the PRD5.5 at time zero after the explosion: 
i.e., before significant radioactive decay has occur. Elements are denoted by solid curves; 
radioactive nuclides by dotted curves. 
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in emission shows solid evidence for the S II "W" (iHamuy et al.ll2003l ). If S II is used as 
the defining feature of SN la, then the synthetic spectrum of Fig. [2] is not that of a SN la. 
Since the angle averaging will automatically lead to some washing out of the flux spectrum 
we should be careful not to draw too strong conclusions about the absence of S II lines, but 
their complete absence in the presence of Si II is interesting. While the relative amount of 
silicon and sulfur in the model is quite similar S III is less abundant than Si III only by 
a factor of two or more everywhere in the outer atmosphere, thus it is unclear why there 
is no evidence for S II in the synthetic spectra (the lines of the ionization stage just below 
the most abundant one tend to be strongest, since there is then energy to excite the levels 
above the ground state). We speculate that the complete absence of evidence for S II in 
the synthetic spectrum is a peculiarity of the particular NLTE state realized in the outer 
ejecta of PRD5.5. If so (and only a very detailed decomposition of the spectrum formation in 
future work can tell), only a PHOENIX level calculation could reveal this peculiarity inherent 
in PRD5.5. 

Figure [3] shows the best fit spectrum compared to SN 2006D seven days prior to B 
maximum, where the C II A6580 line is prominent in the observed spectrum. At this earlier 
time the emission portion of the carbon line in the synthetic spectrum is evident but the 
absorption trough is not blue enough and does not have the right shape. The triangular shape 
in the observed spectrum is almost certainly due to line blending. S II is again missing in 
the synthetic spectrum. The outermost compositions used in the spectral modelling are 
somewhat uncertain due to lower resolution in the outer parts of the spectral grid and to 
the need to interpolate from the gridless hydro models (the grid needed for good resolution 
in synthetic spectra is different from that needed in hydro). This effect, coupled with the 
smearing created by the angle averaging procedure would lead us not to expect perfect 
agreement with the observed lineshapes. 

One of the key features of the pulsating-reverse-detonation model is that it creates a 
composition inversion. That is, the initial deflagration burns at high density and thus tends 
to make iron-peak elements on the outside at high velocity. High velocity features have been 
detected in supernova spectra at early times, but they tend to be either unburn t material 



such as carbon, or intermediate mass elements like calcium. [Branch et al.l (120061 ) used high 
velocity Fe II to fit near maximum light supernova spectra, using the simple parameterized 
code SYNOW. However, the optical depths used were low, and thus the material could be 
primordial. It is impossible to use SYNOW to directly determine abundances. Figure H] 
shows the combined HST and optical spectrum of SN 1992A five days past maximum light 



(IKirshner et al.lll993l ). Plotted on this figure are five different choices for the outer boundary 
condition, the total luminosity in the observer's frame. The luminosity varies from the lowest 
value to the highest by a factor of 13. The logarithmic bolometric luminosity (ergs s~^) is: 
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2000 3000 4000 5000 6000 7000 8000 9000 

A (A) 

Fig. 2. — The synthetic spectra of a full NLTE model of PRD5.5, 20 days after explosion, 
is compared to the observed spectrum of SN 1994D on March 21, 1994 (the time of B 
maximum). The observed spectrum has been corrected for redshift assuming a velocity of 
448 km s~^ and a reddening of E{B — V) = 0.06. 
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2000 3000 4000 5000 6000 7000 8000 9000 

A (A) 

Fig. 3. — The synthetic spectra of a full NLTE model of PRD5.5, 13 days after explosion, 
is compared to the observed spectrum of SN 2006D on 7 days prior to B maximum. The 
observed spectrum has been corrected for redshift z = 0.00853. No reddening correction was 
made. 
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Ll=-17.02, L2=-18.58, L3=-19.07, L4=-19.31, L5=-19.81. Clearly no value chosen for the 
luminosity can reproduce the observed colors across the UV+optical spectrum. This is 
almost certainly due to the large amount of iron-peak elements in the outer parts of the 
model. These elements have such strong line blanketing that they always push much of the 
flux to the red and thus the red colors are too red for all reasonable luminosities. 



5. Conclusions 



The pulsating-rcvcrsc-dctonation model is a novel mechanism that combines the advan- 
tages of the delayed detonation mechanism (higher explosion energy, low amounts of unburnt 
fuel remaining) without requiring a deflagration to detonation transition in an unconflned 
medium. However, the mechanism still has to be explored further and different variations 
will change the total mass in the initial deflagration. We have examined just one particular 
realization here and we have also studied only the angle averaged flux spectra calculated in 
spherical symmetry, thus losing some of the important 3-D effects of the model. We note 
that we have also studied the spectra of PRD6.0. PRD6.0 burned carbon more efficiently 
during the deflagration phase, thus the mass fraction of stable iron and nickel in the outer 
layers of that model is nearly one order of magnitude lower than in PRD5.5. However, 
the mass fraction of ^^Ni in the outermost layer is only slightly lower (~ 15%) in PRD6.0. 
The overall differences in the synthetic spectra were not significant. This suggests that the 
composition inversion inherent in the model will lead to supernova spectra and light curves 
that are extremely red, which has not been observed. A way out of this would be if there 
were some pre-expansion during the initial deflagration phase so that the burning is to in- 
termediate mass elements rather than all the way to the iron-peak. The pre-expansion is 
dependent on details of the initial model. For instance, if a shallow enough thermal gradient 
were present in the bubbles, the nuclear burning would propagate initially with a shockless 
supersonic phase velocity that would lead to pre-expansion before the combustion settled in 
a subsonic flame regime. Another possibility is that the explosion takes place in the core 
of a rotating white dwarf. Since for the same central density, fast rotators are more prone 
to synthesize more intermediate-mass elements and less iron than static initial models. The 
determination of the ignition conflguration of the white dwarf is one of the key problems in 
the physics of SNIa during the explosion phase. Thus, we reiterate that the problem found 
in the models studied here do not rule out the PRD scenario, just the particular realizations 
that we have studied. Clearly there is more work to be done in studying this particular explo- 
sion mechanism. Our work also shows the power of quantitative spectroscopy in validating 
hydrodynamical explosion models. 
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2000 3000 4000 5000 6000 7000 8000 9000 

A (A) 

Fig. 4. — The synthetic spectra of a full NLTE model of PRD5.5, 25 days after explosion, are 
compared to the observed spectrum of SN 1992A observed five days past B maximum. The 
indicated lines show synthetic spectra of models with increasing total bolometric luminosity 
in the observer's frame (L1-L5). The observed spectrum has been corrected for redshift 
assuming a velocity of 1845 km s^^. No reddening correction was made. While the optical 
spectra can be somewhat fit, no variation in the luminosity can fit the entire wavelength 
range from the UV to the optical. 
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